Thermal Canting of Spin-Bond Order 
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Magnetism arising from coupled spin and spatial degrees of freedom underlies the properties of 
a broad array of physical systems. We study here the interplay between correlations in spin and 
space for the quantum compass model in a finite external field, using quantum Monte Carlo meth- 
ods. We find that finite temperatures cant the spin and space (bond) correlations, with increasing 
temperature even reorienting spin correlations between orthogonal spatial directions. We develop 
a coupled mean field theory to understand this effect in terms of the underlying quantum criti- 
cal properties of crossed Ising chains in transverse fields and an effective field that weakens upon 
increasing temperature. Thermal canting offers an experimental signature of spin-bond anisotropy. 

PACS numbers: 75.10.-b, 03.67.Lx 
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INTRODUCTION 

The collective behavior of spins in highly correlated 
quantum magnets typically establish order globally, over 
a large portion of the system by breaking a symmetry 
of the system. Conventional characterizations of mag- 
netic order, such as ferromagnetism, therefore rely on 
microscopic models with global symmetries. Recent the- 
oretical studies have, in contrast, analyzed models that 
support quantum phases without long-range ordering, for 
which novel types of order, including intriguing topologi- 
cal quantum liquid phases, can nevertheless arise [2, 0, y[ . 
Quantum liquids can be difficult to identify (both theo- 
retically and experimentally) because the lack of long- 
range order defies characterization by simple bulk order 
parameters such as magnetization. Models with quantum 
liquid ground states can harbor quasi-local symmetries 
which impose a symmetry on a small, non-extensive sub- 
set of spins Q ■ Examples of such quasi-local operators in- 
clude one-dimensional operator chains embedded in two- 
dimensional graphs. The qualitatively distinct character- 
istics of systems with reduced, non-extensive symmetry 
play an important role in identifying novel quantum liq- 
uid phases. 

We study here one of the simplest of these models, the 
quantum compass model in a magnetic field, h: 
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where Sfj 



erfj/2 are spin 1/2 operators at the [x, z) 
coordinates on an L x L square lattice. It is suffi- 
cient to consider 7 > 0, since with a tt rotation of spins 
about the z-axis on a single sublattice, our results apply 
to the 7 < case as well. The novel coupling between 
bond and spin degrees of freedom separates this class of 
models from more conventional models of magnetism and 
motivates the intriguing question as to whether there arc 
generic signatures of anisotropic spin-bond coupling. 



The quantum compass and related models were first 
discussed in the context of Mott insulators [H| and have 
been studied as simple models of orbital order in the 
transition metal oxide compounds with anisotropic cou- 
pling among orbitals defining pseudo-spins @, @, 0] • For 
example, recent experiments [3 0, E3] on a two dimen- 
sional e g orbital compound, LaSrMn04, observe intrigu- 
ing thermal effects. Anomalies in thermal expansion 
measurements Q and in Raman scattering [Toj . as well 
as structural changes measured by neutron and x-ray 
diffraction have all been interpreted as indicating a 
change of orbital occupation that is driven by tempera- 
ture. The reorientation of orbital direction with increas- 
ing temperature implied by the experiments raises the 
question as to whether this results from thermal repopu- 
lation of a non-interacting orbital system, or from com- 
petition between orbital-orbital interactions and thermal 
effects. Detailed understanding of this phenomenon re- 
quires quantitative analysis of the competing roles of su- 
perexchange, phonon, crystal field and Jahn- Teller ef- 
fects. However, the orbital-only degrees of freedom in 
these systems have been shown to be approximately de- 
scribed by a quantum compass model with anisotropic 
pseudo-spin interaction terms given by superexchangc 
contributions and an effective magnetic field determined 
by the crystal field splitting . This enables one to ad- 
dress in general the intriguing question of whether and 
how orbital-orbital interactions can influence the ther- 
mal redistribution of orbitals. We answer this question 
in the affirmative here, showing that a novel type of ther- 
mal reorientation, 'thermal canting', of nearly degenerate 
orbitals can indeed arise from a large anisotropic orbital- 
orbital interaction of the type exemplified in Eq. [TJ We 
also show that this thermal canting can provide an ob- 
servable signature of the anisotropic interactions. 

Eq. ([T]) has been studied in other contexts as well. A 
duality mapping II 21] relates Eq. ([T]) to models of p-w&ve 
superconductors j 1 31 ] - Recent proposals put forward to 
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realize Eq. |T]) in optical lattices, using asymmetric tun- 
neling |l4j or polar molecules flEl ]. and in Josephson junc- 
tion arrays [l6[ (finite sized versions of which have re- 
cently been realized have been motivated by the 
recognition that the h = limit is characterized by sets of 
low energy two-fold degenerate states that may provide 
protected subspaces for encoding quantum information 
ljj |. In addition, related [H anisotropic spin mod- 
els (Eq. (4) in 0]) show topological order which allow a 
robust form of quantum information processing 0, [H| . 
A generic experimental signature of the spin-bond asym- 
metric models would be a valuable tool in the search for 
realizations of models with anisotropy and related quasi- 
local symmetries, a In this paper we study the finite 
temperature properties of the quantum compass model 
for h > 0, using both quantum Monte Carlo (QMC) 
and mean field theory to identify generic signatures of 
the novel set of one-dimensional symmetries underlying 
Eq. (p}. We find a remarkable thermal canting effect 
whereby finite temperatures enhance spin-bond correla- 
tions r x , where r a = L~ 2 J2ij{S"jS" +1 j) and a = x, z. 
The enhancement resembles an order-by-disorder process 
@j 0, [2l|. The reorientation of net spin-bond correla- 
tions (from r z to r x ) defines a phase akin to a liquid 
crystal smectic C phase where molecular moments order 
along a vector ('director') tilted with respect to transla- 
tionally ordered columns. We develop a mean field theory 
using intersecting Ising models and their crossover phase 
diagrams [22| to show that the finite temperature reduc- 
tion of the z-magnetization leads to thermal canting at 
the mean field level and should therefore play a role in 
similar models. For 7 > 1 we find that finite tempera- 
tures can completely reorient spin-bond correlations to 
point in orthogonal directions in spin and real space. 



SYMMETRIES IN ZERO FIELD LIMIT 

We begin our study with an analysis of the unique sym- 
metries of the quantum compass model in several well- 
defined limits, for h = 0. Under high strain, 7 3> 1 
(7 <C 1), the system forms a series of L weakly cou- 
pled spin chains lying along the x (z) direction in real 
space. At zero temperature the chains order to point 
along the x (z) directions in spin space. Note that, with 
h = 0, Eq. |T]) commutes with the quasi- local symme- 
try operators Pi = TTj %Sf,j and Q 3 = fJi > which 
act on ordered chains to generate spin flips along z and 
x chains, respectively (P is not an exact symmetry for 
h 9^ 0). Gapped excitations along chains are then finite 
length domains of flipped spins, as in one-dimensional 
Ising models. The operators Q (P) can flip all spins along 
the ordered chains in the ground state, resulting in yet 
another ground state. This results in a two-fold degen- 
eracy of each chain and therefore a massive ground state 
degeneracy of the two-dimensional system of at least 2 L . 



As 7 — > 1, the inter-chain coupling allows a set of inter- 
chain excitations. Ref. [23| showed that high order inter- 
chain fluctuations of the magnetization cost zero energy 
(for L — > 00), preserving a large, 0(2 L ), ground state 
degeneracy for all 7. The 7 = 1, h = limit yields 
discrete, global rotational symmetries (see Ref. [12| for 
a discussion). A nematic-like order parameter has been 
defined 13, 3]: T = r z — r x , to characterize smectic- like 
bond ordering along the x and z directions in spin and 
real space, analogous to the directional ordering of molec- 
ular moments along columns in liquid crystals, r = 
implies a disordered phase with no preferred direction. 
Finite temperature studies 0, HH suggest a phase transi- 
tion from an ordered, \r\ > 0, phase to a disordered phase 
with T c ps 0.055 [H] in units of the interaction strength 
Eq. Example interaction strengths in correspond- 

ing experimental systems are ~ lOOOif for orbitals [H or 



IK for Josephson junctions [17 1 



QUANTUM MONTE CARLO SIMULATIONS AT 
FINITE FIELDS 



A weak finite magnetic field picks out a single ground 
state with a gap to excitations set by h, simplifying QMC 
simulations. We employ the stochastic series expansion 
using directed loops [26 



code in the ALPS 




A modification of the 



package allows a treatment of 
bond terms in the updating scheme. All QMC 
results presented here are converged to the thermody- 
namic limit in studies of several system sizes up to L = 30 
for periodic boundary conditions. 

Our results show that competition between the inter- 
action and the magnetic field sets up an interesting in- 
terplay in the energetics along both spatial directions. 
The insert in Fig. Q] (top panel) shows a schematic of 
the generic behavior observed in the QMC calculations 
for 7 < 1 as we increase temperature. At low tem- 
peratures the weak magnetic field aligns the net mag- 
netization along the z-direction. The magnetic field 
and z-interaction term favor spin chains that are bond- 
correlated along the z direction, i.e. large r z . Using mean 
field theory we will show that inter-chain interactions 
contribute to an effective magnetic field that suppress 
excitations along the x-dircction. As we increase tem- 
perature above the gap set by h, the applied magnetic 
field becomes less relevant while the mean field effective 
magnetic field is reduced with increasing temperature. 
Excitations along x are then allowed to recover. Conse- 
quently, bond correlations build up along the x-direction 
yielding r x — > r z as can be seen in the QMC and mean 
field results in Fig. Q] At high temperatures both types 
of correlations are trivially suppressed as the system be- 
comes thermally disordered. 
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Figure 1: The top (bottom) panels plot the spin-bond correla- 
tions versus temperature along the x (z) direction computed 
via quantum Monte Carlo (open circles) and partial mean 
field theory (solid line) for the quantum compass model on a 
square lattice in a magnetic field. The peak in r x indicates 
a reorientation of bond correlations with increasing temper- 
ature depicted schematically in the insert, where j and — > 
correspond to spins aligned along the z and x directions, re- 
spectively. The data are plotted here for 7 = 1 and h = 0.08. 



PARTIAL MEAN FIELD THEORY 

At first sight, the build up of bond correlations along 
the x-dircction is puzzling since one would expect a sup- 
pression of bond correlations with increasing tempera- 
ture. To understand this build up of correlations along 
the x-direction, we construct a coupled partial mean 
field theory (PMFT) in which spin interactions along a 
given bond direction are treated at the mean field level 
while the interactions along the orthogonal direction are 
treated exactly. In our two-step protocol we first choose 
the z direction (indicated by (...) a ) and make a mean field 
decoupling, S*jS?j +1 — > S* t j{S*j +1 ) z to generate the ef- 
fective field Bf = {S z J+1 ) z + (S z , j _ 1 ) z + h for spins along 
x. We then repeat the process by exchanging x and z, 
resulting in a set of two coupled mean field equations: 

H MF = - X] S tj S i,j+l ~ k ^2 S t] ~ B * S ^2 S i,ffl 

h mf = -I'yi^ti^i+i.i ~ B t s y^; s ti ( 3 ) 

While in general one may take Bf = -y{(Sf +hj ) x + 
(Sf_i j) x ) and make a self-consistent solution to Eqs. ([3]) 




Figure 2: Crossover diagram for a one-dimensional transverse 
field Ising model. The ic-axis indicates the effective magnetic 
field which can be lowered to tune from the paramagnetic 
side (weak nearest neighbor correlations), through a quantum 
critical point, B c = 7/2, to a regime with relatively strong 
nearest neighbor correlations. The dashed lines schematically 
define the quantum critical regime. For the PMFT defined 
by Eq. ([3J, Bf is temperature dependent and plotted for 
h = 0.08 (solid line) and h = 0.001 (dotted line) with 7=1. 
Inset: The same as Fig.[T]but for magnetization density. The 
PMFT curve plots (S*) x . 



and , in what follows we keep the mean field solutions 
analytic by setting Bf = 0. This approximation does 
not qualitatively alter our results for 7 near unity. We 
also assume a uniform net magnetization (S z ) z = {Sf j) z . 
The magnetization (S z ) z obtained from solving Eq. © 
can then be substituted directly into Eq. ([3]). This mean 
field approach is motivated by the work of Ref. [3^] but 
differs by i) having a set of two coupled mean field equa- 
tions, ii) addressing a new parameter regime of finite tem- 
peratures and magnetic fields, and iii) obtaining qualita- 
tively accurate results without a self consistency loop. 
The mean field approach is expected to be qualitatively 
accurate because the finite magnetic field suppresses fluc- 
tuations along the z direction. 

We can now solve the PMFT equations analytically. 
Solving Eq. ^ using the well known solutions of the clas- 
sical Ising model |3l(, we obtain (S z ) z = tanh(C)/(2G), 
where G = [1 + (exp(-4RT) - l^ech^C)] 1 / 2 , K = 1/4T, 
and C = h/2T and can estimate r z as: (Sf jSf j + \) z = 
(l/4)d A -ln{exp(f0cosh(C)(l + G)}. This first estimate 
for the magnetization is then used as input to Eq. ([3J, 
which is the Ising chain in a transverse field and can 
be solved exactly via the Jordan- Wigner transforma- 
tion [33]. We obtain the total energy: {Hf [F ) x = 
Ek^(nk ~ 1/2), where e k = [(Bf) 2 - jBf cos (fc) + 
7 2 /4]V 2 anc j nk i s the Fermi distribution function. 
Our PMFT estimate for the magnetization of the two- 
dimensional system becomes: {S z ) x = J^kftBf — 
7COS (fc)](4efc) _1 tanh (e^/2T). We can also obtain an 
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Figure 3: Net spin-bond correlations as a function of temper- 
ature for h = 0.1 and 7 = 1/0.85. Thermal canting reorients 
the spin-bond correlations to lie predominantly along the x 
direction in both spin and space for T > 0.25. 



estimate for r x : (SfjSf +1J ) x = -(l/jL){Hl IF ) x - 
(l/"/)Bf{S z ) x . By setting Bf = we ignore the in- 
fluence of the magnetization along x on r z . The result is 
an over-estimate for r z at low temperatures. Fig. [I] and 
the inset of Fig.[J]show a comparison between the PMFT 
and QMC results for the temperature dependence of r x , 
r z , and (S z ). It is evident that the PMFT captures the 
qualitative features of the QMC results. 

The success of the PMFT analysis implies that finite 
temperature excitations of the two-dimensional quantum 
compass model can be approximated by excitations along 
crossed (coupled) one-dimensional Ising models. The 
peculiar enhancement of r x with temperature can now 
be understood from the crossover phase diagram of the 
transverse field Ising model, Eq. ([3]) . Figure [5] plots a 
schematic of this crossover phase diagram. In our mean 
field analysis, B x carries a temperature dependence be- 
cause the z-spins demagnetize with increasing tempera- 
ture. In Fig. [2] we superpose B X S (T) from a solution of 
the mean field equations, which shows that an increase 
in temperature lowers the effective magnetic field seen by 
spins along the x direction. The bond correlations along 
the x direction are therefore enhanced because of a con- 
sequent drastic reduction of {S z ) z . We conclude that the 
nonlinear temperature dependence in the effective mag- 
netic field experienced by the spin chains lying along the 
x direction is therefore responsible for the thermal cant- 
ing of spin-bond correlations observed in QMC. 



ft~7>l limit of Eq. (p} can be thought of as a set of 
nearly independent transverse field Ising models with a 
quantum critical point (B c — 7/2). Lowering h takes the 
system from the paramagnetic phase through the quan- 
tum critical point (at T = 0) to an ordered chain phase 
along the x direction. For a specific choice of 7 > 1 
we can therefore completely reorient bond correlations 
by tuning h. Thermal canting also tunes B x but with 
a directionally independent parameter, temperature. At 
low temperatures the finite magnetic field favors bond 
correlations along the z-direction, r > 0, while at higher 
temperatures thermal canting of spin-bond correlations 
now favors a net orientation along the x-direction, r < 0. 
This implies a canting transition at some temperature 
Tcant, which is illustrated in Fig.[3]for 7 = 1/0.85, where 
Tcant ~ 0.25. These results predict a striking effect, 
namely, a complete reorientation of net bond correlations 
in both spin and real space with increasing temperature. 
This thermal reorientation arises from a distinct mecha- 
nism, thermal canting, that can again be understood in 
the PMFT calculations as being due to a weakening of 
the effective magnetic field with increasing temperature: 
thus it is a consequence of the two-dimensional, vector 
nature of the spin-bond correlations in Eq. ([1]) . 

Our two-dimensional spin-bond analysis at small finite 
h also provides a microscopic rationale for the |r I — 
transition observed in the h = QMC simulations 0, [SH 
in terms of an immediate thermal canting of spin-bond 
correlations. Our QMC and PMFT results both show 
that the r x peak in Fig. [1] moves to lower temperatures 
with decreasing magnetic field. The mean field analysis 
indicates that this enhancement of r x at low tempera- 
tures derives from the rapid depolarization of (S z ) z with 
increasing temperature when h is small. Mean field re- 
sults for extremely small h values and 7 = 1 are shown in 
Fig. [4] We see that the increase in r x with temperature at 
low h concomitantly tends to decrease r (see inset). Ex- 
trapolating the PMFT results to h — ► yields, for 7=1, 
r x = (l/4)tanh(l/4T) + e>(fr 2 ) and r = + O{h 2 ). Thus 
r disappears altogether for h — > 0. This contrasts with 
the h = QMC results of Ref. @, HI] where a low tem- 
perature plateau in |r| was followed by a rapid decrease 
at larger T and assigned to a disordering transition. The 
difference may be due to a lack of fluctuations in the 
PMFT which could play a large role in the h — > limit. 



SUMMARY 



MAGNETIC FIELD DEPENDENCE OF 
THERMAL CANTING 

We now examine thermal canting in finite fields h in 
greater detail. Fig. [3] indicates that the net spin-bond 
correlations completely reorient with increasing temper- 
ature. To understand this behavior, first note that the 



We have found a distinct, thermal canting effect in 
the x-z plane of the two-dimensional quantum compass 
model. Excitations governed by the rather unique set 
of chain symmetries implicit in the model provide a mi- 
croscopic mechanism responsible for this effect. Thermal 
excitations above the paramagnetic gap can be approxi- 
mated by those of intersecting Ising chains that are gen- 
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Figure 4: PMFT results for the bond correlations along the x 
(solid line) and z (dashed line) directions versus temperature 
for 7 = 1 and several values of h. The inset shows the behavior 
of the scalar difference r = r z — r x . 



crated by quasi-local operators stemming from the exact 
symmetry, Q, and the approximate symmetry, P. Such 
excitations lead to an enhancement of spin-bond correla- 
tions along the x direction at higher temperatures. This 
study shows that temperature can be used as a parameter 
to experimentally tune the balance between the excita- 
tions generated by these operators and to thereby reveal 
the underlying chain symmetries and anisotropic inter- 
actions. Thermal canting arises in our mean field theory 
(as well as QMC). Our partial mean field argument ap- 
plies to a number of other lattice geometries, in addition 
to the rectangular lattice of Eq. [TJ 

In summary, our combined QMC and partial mean 
field analysis show that anisotropic interactions can lead 
to a thermal redistribution of spin-bond correlations. 
The results suggest that observation of thermal cant- 
ing can provide an indicator of spin-bond anisotropy, 
with significant implications for the transition metal ox- 
ide compounds that are described by pseudo-spin mod- 
els of orbital degrees of freedom. Our demonstration that 
thermal canting can be driven by anisotropic interactions 
suggests that the recent observations of a thermal re- 
orientation of nearly degenerate orbitals in LaSrMn04 
BI0 warrant further analysis in terms of anisotropic 
pscudospin models of orbitals. 
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